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The  effort  as  summarized  in  the  report  provides  accurate  theoretical  analyses 
of  the  optical  power  transmission  properties  of  a  number  of  devices  important 
to  RADC  efforts  under  TPO  H(D  -  Solid  State  Devices.  Subthrust  #3  -  Electro- 
Optical  Components.  The  devices  include  couplers,  tapers,  horns,  and 
branches.  The  programs  anu  techniques  provided  by  ttte  contractor  permit 
in  effect  computer  experiments  to  be  done  for  a  very  large  variety  ot  design 
parameters.  To  do  the  actual  experiments  in  the  laboratory  with  this  range 
of  parameters  would  be  enormously  more  expensive  and  time  consuming. 

The  report  is  then  a  crucial  element  in  simplifying  and  accelerating  the 
design  process  and  in  leading  to  final  design  specs  in  the  shortest  time. 


Project  Engineer 


I .  INTRODUCTION 

This  final  report  summarizes  the  work  performed  under  Contract  F19628-80- 
C-0053  which  the  Electronic  Systems  Division  of  the  Air  Force  Systems  Command 
granted  to  the  EMtec  Engineering,  Los  Angeles,  California.  The  work 
was  begun  in  January,  1980  and  completed  in  August,  1980. 

The  principal  thrusts  of  this  R  &  D  study  in  performing  numerical  analysis 
of  multimode  fiber  components  were  two  fold.  Firstly,  we  wish  to  learn  the 
limitation  (and  possible  improvement)  of  our  numerical  scheme^and  secondly, 
we  wish  to  obtain  numerical  data  for  realistic  multimode  fiber  structures. 
Specifically,  the  following  tasks  were  carried  out: 

a)  Study  the  effect  of  step  index  gradient  and  of  tight 
beam  confinement  by  an  adaptive  coordinate  scheme. 

b)  Study  the  effect  of  the  presence  of  absorber  at  the 
edge  of  the  mesh  on  the  beam  propagation  character¬ 
istics  of  multimode  fiber  structures. 

c)  Compute  the  coupling  characteristics  of  tapered  multimode 
fiber  couplers  and  unequal  size  fiber  couplers  . 

d)  Obtain  data  for  reflection  coefficients  and  beam  waist 
changes  for  multimode  fiber  tapers,  horns  and  branches. 

In  section  II  we  shall  present  the  implementation  of  the  adaptive 
coordinates  in  our  numerical  solution  of  the  scalar  wave  equation.  Then, 
the  scheme  to  include  an  absorber  at  the  edge  of  the  mesh  will  be  de¬ 
scribed.  Finally,  an  approximate  approach  to  obtain  the  reflection 
coefficients  for  complex  fiber  structures  will  be  shown.  Detailed  re¬ 
sults  of  our  study  on  the  proposed  tasks  are  given  in  Section  III. 

Concluding  remarks  and  recommendations  for  future  work  are  included  in 
Section  IV. 
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II.  ANALYTICAL  APPROACH 

The  basic  approach  taken  to  find  the  solution  of  wave  propagation  along 

complex  f^ber  structures  is  to  solve  the  reduced  scalar  wave  equation  via  the 

2 

fast  Fourier  transform  (FFT)  technique.  In  this  section  we  shall  first  indicate 
the  conditions  underwhich  the  exact  vector  wave  equation  may  be  simplified  to 
yield  the  reduced  scalar  wave  equation.  Then  we  shall  introduce  the  concept 
of  adaptive  coordinates^  and  incorporate  this  concept  in  the  solution  of  the 
reduced  scalar  wave  equation  via  the  fast  Fourier  transform  technique. 

A.  Formulation  of  the  Scalar  Wave  Approach.  Starting  with  the  vector 
wave  equation  for  the  electric  field  vector  E  in  the  fiber  structure, 

2 

VxVxE-w^E-O  (1) 

where  u)  is  the  frequency  of  the  wave,  the  permeability  and  and  c  =  f(r),  the 
inhomogeneous  permittivity  of  the  structure,  and  making  use  of  the  vector 
identity 


V  x  V  x  E  =  V(V»E)  -  V  E 


(2) 


and  the  relation 


V* E  =  -  , 


(3) 


one  has 


V2E  +  W2  (1Q( E  -  P(yVf-E)  =  0 

Rewriting  Eq.  (A)  gives 


(A) 


2  2 
V  E  +  oj  ft 


0<0  f.'0 


~~  V  (  7  V'*E> 

o 
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The  relative  importance  of  the  terms  within  the  curly  brackets  can  be 
determined  from  the  following 


where  the  symbol  0  means  the  "order  of  magnitude,”  and  /  is  the  smaller  of 
the  distance  over  which  */(q  and  E.  change  appreciably.  For  single-mode 
fiber  structures,  the  values  of  ( /(q  and  kQA  are  typically  in  the  range 

</(Q  «  *(2)  (8) 

k Ql  1=  f  (1C)2  or  1q3)  >  I*  to  100  )  (9) 

A  =  0(lfi) 

It  follows  that  the  second  term  within  the  curly  brackets  in  Eq.  (5)  is 
several  orders  of  magnitude  smaller  than  the  first  term  (/(q  E.  It  is  there¬ 
fore  justifiable  to  neglect  the  second  term  and  write  Eq.  (5)  in  the  form 

2  2  e 

v  E  +  k  —  E  “  0 
—  o  e  — 
o 

The  physical  significance  of  replacing  Eq.  (5)  by  Eq.  (10)  is  this.  By 
discarding  the  term  ? -  V?  E,  we  are  neglecting  any  depolarization  effects 
that  may  occur.  This  means  that  the  wave  retains  the  polarization  it  has 
at  the  source,  which  is  evidenced  by  the  fact  that  Eq.  (10)  can  be  reduced  to 
a  scalar  equation  by  writing  E(x)  in  the  form 
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h(x)  “  ®pU(x) 


where  e  is  a  unit  vector  in  the  direction  of  the  initial  polarization  of  the 
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wave.  Subsituting  Eq.  (11)  in  Eq.  (10),  we  find  that  u(x)  satisiies  the  scalar 


wave  equation. 


7  2  r 

V  u  +  k  —  u  -  0 

°  £o 


This  equation  with  the  boundary  condition  on  the  initial  surface,  and  the 

radiation  condition  at  infinity,  completely  specifies  u (x) ,  from  which  we  can 

then  obtain  the  electromagnetic  field  vectors  E  and  H. 

ikn  z 

If  we  write  u  as  the  product  of  a  factor  e  that  accounts  for  the 

rapid  change  in  the  phase  of  u  along  the  direction  of  propagation  and  a  complex 
amplitude  A(x,z),  a  further  simplification  of  the  problem  results 


* 


t2ikn  r —  +  V~  +  k^(n2(x,z)  -  n2)  ]  A(x,z)  =  -  a2A(x,z) 


2  c»  3  ^ 

where V  is  the  transverse  Laplacian  — _  +  — r— ,  and  n  is  a  given  constant 

A  n ..  A  Z  U 

3x  3y 

which  represents  the  refractive  index  of  some  uniform  medium.  At  laser  wave¬ 
lengths  the  complex  amplitude  A(x)  varies  much  more  rapidly  transverse  to  the 
direction  of  propagation  than  it  does  along  the  direction  of  propagation.  This 
enables  us  to  make  the  paraxial  approximation  wherein  the  term  on  the  right  side 
of  Eq.  (13)  is  neglected  (  in  the  Russian  literature  this  is  called  the  parabolic 
approximation).  So,  the  complex  amplitude  now  satisfies 


[l2kno  3z 


2  2  2  I 

+  — r  +  k  (n  (x,z)  -  n  A(x,i 

3v2  "  °J  " 
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For  given  initial  data,  i.e.,  values  of  the  field  at  points  on  the  initial 
surface,  the  propagation  simulator  must  generate  the  corresponding  field  values 
at  the  terminal  aperture  such  that  Eq .  (14)  is  satisfied.  To  do  this  we  divide 
the  medium  into  slabs  defined  by  planes  on  which  z  is  constant.  In  going  from 
one  slab  to  the  next,  we  write  A(x,z)  in  the  form 

A(x,z)  =  er^X,Z^  w(x,z)  q5\ 


where  F  (x,z)  is  a  phase  function  associated  with  the  medium  inhomogeneities 


z 

r(x,z)  »  J  £n2(x,y,z')  -n2Jdz'. 


(16) 


The  modified  complex  amplitude  w(x,z)  then  satisfies  the  equation 


[i2kn0  h  +  e~r  eF  ]  w(^’z)  =  0 


(17) 


with  the  initial  condition 

w(x,y,0)  =  u(x ,y,0) 

Physically,  these  equations  approximate  the  propagation  in  the  inhomogeneous 
medium  by  a  two-step  process  at  each  z  increment.  First,  we  propagate  the 
field  u(jc  )  at  z  -  Az/2  to  z  +Az/2,  assuming  that  the  intervening  space  is 
homogeneous.  The  effect  of  the  inhomogeneitics  between  z  -  Az/2  and  z  +  Az/2 
is  then  accounted  for  by  multiplying  this  solution  by  the  phase  factor  exp(r). 


B.  Adaptive  Coordinates 

To  reduce  the  size  of  the  mesh  required  to  solve  Eq.  (17)  numerically. 


61/2  -  a 


*p: 


(23) 


where  is  a  characteristic  dimension  of  the  beam  at  the  initial  surface 
(e.g.,  the  e-fording  radius  of  a  gaussian  beam),  f  is  the  distance  to  the  focus, 
and  a  is  a  constant  determined  by  the  requirement  that  the  solution  be  confined 
within  the  boundaries  of  the  mesh  at  the  focal  plane.  The  choice  a  =  1  yields  a 
coordinate  system  that  converges  at  a  rate  determined  by  the  free-space  dif¬ 
fraction  of  a  gaussian  beam  having  an  e-folding  radius  P^. 

When  written  in  terms  of  the  converging  coordinate  variables  defined  above, 
Eqs.  (15)  and  (17)  for  the  complex  amplitude  are  replaced  by  the  relations 

w(x,y,z)  =  £)(£_,£)  expfr)  v(£,s;  (2A) 


w  (£_,£)  =  (a1^2  N  (z)  j  1  exp  [  J  (cj  +  c,\) 

~r‘%fZZ + 


“  2  exP(-f )  (  2  +  3  2  )  exp(r) 

®  C 1  3^2 


v  =  0 


(25) 


(26) 


(27) 


whereA£  Is  the  increment  in  £  in  going  from  z-Az/2  to  z+A  z/2.  The  initial 
condition  for  v  is 


v(£_ ,£)  =  w(x,y,z)/5  (r,C) 


(28) 


To  solve  Eq.  (27)  we  utilize  the  fact  that  for  sufficiently  small  values 
of  AC  (i.e.,Az)  the  effect  of  the  exponential  factors  exp  (-  r  )  in  this  equation 
is  small.  Hence,  we  solve  the  simpler  equation  obtained  when  these  factors  are 
equated  to  unity 


JL 


0 


(29) 


We  use  a  fast  Fourier  transform  technique  to  solve  Eq.(29).  The  basis 
of  this  approach  is  the  fact  that  the  solution  of  Eq.  (29)  can  be  expressed  in 
the  form  of  a  discrete  Fourier  series 

N-l  N-l 


2  I 

n=r  n  Y\~t 


l'tnn(C'tl  «‘p[‘<PnCi+qnC2>l 


(30) 


where  the  Fourier  coefficients  V  are  determined  from  the  initial  data  and 

mn 

Eq.  (29)  as  follows.  The  initial  values  of  V  are  obtained  by  taking  the  discrete 

mn 


Fourier  transform  of  the  initial  values  of  v(  [  »  C  over  a  mesh  of  points 
Cx  -  U-(-/2)]  A?,  C2  =  l  j“  (N/2)  ]  AC  U.j  =  0,1,...,N-1) 


V  (5;) 

mn  i 


_  (-D 


m+n 


N-l  N-l 


N* 


*-U  j-0 


expj-  1—  (m£  +  nj)| 


(31) 


The  dependence  of  Is  then  determined  by  subsituting  Eq.  (30)  in  Eq.  (31), 


which  yields 


dV  ■  , 

. _ 5HL  +  1  ( 

dZ  2  \ 


p  +  q*  V  =0 
m  /  ran 


from  which  it  follows  that 


V  (f)  =  V  (£  .) 

mnVj'  mn  s 


-  lM  *  qn  )  45 


Finally,  it  can  be  shown  that  in  order  for  the  discrete  Fourier  series  representation 

of  v  given  in  Eq.  (30)  to  be  real  when  v  is  real,  the  coefficients  p  and  q  must 

m  ti 

have  the  form 


„  -  2n  / 

Pm  "  NAC  \  2/ 

2tt  (  N\ 

qn  NAC  \n  _  2/ 


Hence,  for  discrete  points  N/2  )b(, ,  Cg  *  (j  -  N/2)&^(^,  j  *  0,1,...,  N-l) 


{(*-!)  “•(*-!)«•  t) 


=  (- 1)£  +  3 


N-l  N-l 


m=0  n=0 


-  ‘  5  ((~?)  *  ♦  i  £  C  ♦  in 
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I"" W  Iff •  (Ml*.  «P 
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^  2  2  .  » 
where  8  =  2tt  A?/(A£)  •  Note  that  v  is  simply  (-1)  times  the  discrete  Fourier 


transform  of  (-1) 


m+n 


V  <£) 
mn^' 


2  2 

The  effect  of  the  medium  and  the  facto*.  exp[(-i/2fe^  +  T,  ^)A£  ]  introduced  by  the 
coordinate  transformation  is  taken  into  account  at  each  £  step  in  the  calculation 
by  multiplying  the  value  of  v  obtained  in  the  previous  step  by  the  quantity  exp(?*) 
defined  in  Eq.  (26),  i.e.,  the  initial  value  inserted  in  Eq.  (31)  is  exp(T)  times 
the  value  of  v  determined  from  the  previous  steps. 

Using  this  adaptive  coordinate  algorithm  we  have  been  successful  in  our 
treatment  of  various  realistic  multimode  fiber  structures.  Results  are  summarized 
in  Section  III. 

C.  Implementation  of  Lossy  Outer  Boundary. 

It  is  believed  that  the  field  touching  the  outer  boundary  of  the  cladding  region 
of  a  fiber  structure  will  be  attenuated  due  to  radiation  or  absorption.  To 
accomodate  this  situation  in  order  to  further  improve  our  computer  simulation,  we 
have  incorporated  the  presence  of  a  lossy  dielectric  layer  outside  the  cladding 
region  in  our  computer  program.  An  example  of  the  index  profile  of  a  fiber  is  shown 
in  Fig.  1: 


0  <  r  <  a 


a  <  r  <  r. 


n  (1  -  6)  , 
a 


25  Vim 
•  62,5  Pro 


ro<  r 


n  +  in. 
r  i 


Figure  1:  A  Typical  Index  Profile  with  Lossy  Outer  Layer 
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Topical  intensity  patterns  of  beams  propagating  in  a  fiber  with  the  index 
profile  given  by  Fig.  1  are  shown  in  Fig.  2.  It  was  found  that  the  propagation 
characteristics  of  the  guided  beams  are  not  significantly  affected  by  the 
presence  of  a  lossy  outer  layer,  except  when  the  spot  size  of  the  beam  is  larger 
than  the  core  diameter,  as  expected. 


D.  A  Heuristic  Approach  in  Obtaining  the  Reflection  Coefficient. 

One  of  the  a'prior  assumption  in  the  development  of  the  FFT  scalar  wave 
approach  is  that  only  paraxial  rays  are  allowed  and  no  reflection  is  permitted. 

This  assumption  enables  us  to  develop  an  algorithm,  there- by,  we  may  obtain  the 
propagating  field  by  a  forward  stepping  process  as  described  earlier.  As  the  field 
evolves  from  one  z  plane  to  the  next  z  +Az  plane,  the  averaged  value  of  the 
refractive  index  as  seen  by  the  field  may  be  different  as  illustrated  in  Fig.  3: 

.t  ,,t  A 

f* T-  7  +  -a 


A  vs-*,  a  a  c \  v. 

!  i  ■*  I 


iJ  i 

I-  \-‘1  V  <. 


'  .A  ■ 


=  ^5 


Figure  3:  Illustration  of  the  averaged  n  as  seer,  by  the  beam. 


In  affect  one  may  postulate  that  the  wave  is  experiencing  reflection  in  a  medium 
with  longitudinally  slowly  varying  refractive  index  as  shown  below: 


reflected 

wave 


inc ident 
wave 


n(z) 


z 


transmitted 

wave 


n(z)  is  given  by  Fig.  3. 
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Figure  4:  Equivalent  Index  Profile 


The  reflection  coefficient  for  a  plane  wave  propagating  in  this  longitudinally  non- 

uniform  medium  may  be  obtained  according  to  a  formula  derived  for  the  case  of  plane 

6 

wave  propagation  in  stratified  layered  medium: 

r  (z)  =  -  exp  [ -is(z) ]  J  y(z)  exp  lis(z)l  dz 


s  ( z ) 


J 


B(z)dz 


8(z)  =  kQn(z) 


Y  (z)  -  /  28 

n 


This  is  the  heuristic  approach  that  wc  shall  use  to  calculate  the  reflection 
coefficient  for  waves  in  our  multimode  fiber  structures. 
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III.  Results 

The  algorithms  detailed  above  have  been  implemented  in  our  computer  programs. 
Results  for  the  proposed  tasks  are  given  in  the  following: 

(a)  Effects  of  Step  Index  Gradient  on  the  Propagation  Characteristics. 

The  purpose  of  this  study  is  to  learn  the  effects  of  step  index  gradient 
on  the  propatation  characteristics  of  waves  in  a  multimode  fiber  guide.  Let 
us  introduce  the  following  index  profile: 

/  \  t  /r\2m  (o<r<a) 

0  a 

n(r)  =  nQ  -  <5  (a<r) 

where  n^,  m  and  <5  are  given  constants  and  a  is  the  core  radius  of  the 
fiber.  For  a  typical  parabolic  index  profile  fiber,  one  has 

nQ  =  1.48,  6  *  0.02,  a  =  25pm,  m  *  1. 

The  constant  6  must  necessarily  be  small  so  that  the  depolarization  effect 
may  be  ignored  and  the  scalar  wave  approach  may  be  justified.  By  varying 
m,  the  steepness  of  the  index  gradient  may  be  varied  as  shown  in  Fig.  5. 

It  should  be  kept  in  mind  that  even  when  the  FFT  technique  is  capable 
of  handling  steep  index  variations,  the  slope  of  the  index  profile  must 
still  be  gentle  enough  so  that  the  gradient  term  in  the  exact  wave 

equation  (Eq.  (5))  may  be  ignored.  We  have  carried  out  propagation  cal¬ 

culation  for  the  following  specific  cases:  n^  =  1.48,  6  =  0.02, 
a  ■  25pm,  m  *  1,4,6,10.  Higher  m  values  means  steeper  index  gradient. 

As  shown  in  Figs.  6,  no  computational  difficulties  were  encountered  for 

even  the  steepest  case  (m  *  10)  in  which  the  index  changes  from  n  »  1.48 

to  n  •  1.46  in  3pm  distance  for  50  pm  core  diameter  fiber.  However,  one 
should  be  aware  that  we  are  pushing  the  limit  of  validity  for  the  scalar 
wave  approach. 
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From  the  results,  it  is  of  interest  to  note  that  as  m  increases  from  1,  i.e, 
as  the  index  profile  deviates  from  the  parabolic  profile,  the  beam  profiles  no 
longer  remain  to  be  of  gaussian  shapes,  but  take  on  ring-type  structures.  This 
implies  that  the  phase  front  of  the  multimode  beam  is  no  longer  a  monotonic 
function  of  the  radial  distance  but  has  become  an  oscillatory  one. 

We  may  conclude  from  these  calculations  that  our  program  is  capable 
of  handling  problems  with  steep  index  gradient.  The  index  transition 
may  occur  in  a  distance  as  small  as  4X  where  X  is  the  free-space  wavelength. 

The  limiting  factor  apparently  is  the  justification  for  the  elimination  of 
the  depolarization  term  in  Eq.  (5). 

(b)  Beam  Propagation  in  a  Ring  Fiber 

A  typical  single-mode  fiber  has  a  core  diameter  of  the  order  of  10pm. 
Consequently  it  is  more  difficult  to  handle  than  a  multi-mode  fiber.  An 
idea  to  enlarge  the  single-mode  fiber  has  recently  been  put  forth  by  Dr. 

L.  Eyges  of  RADC.  He  suggested  that  perhaps  a  ring-type  structure  may 
support  a  single  mode  and  yet  possesses  larger  dimension  than  the  usual 
solid-core  fiber.  This  task  was  undertaken  to  investigate  this  possibility. 

Let  us  postulate  that  the  index  profile  of  a  ring  fiber  takes  the  following 

^onn  n  =  n  -  6  for  a-A<r<a+A 

0  A 

n  =  nQ  -  <5  for  a-A>r,  r>a+A 


where  n^,  6,  a,  A,  and  m  are  given  constants.  By  increasing  the  m  value, 
one  may  adjust  the  steepness  of  the  index  gradient  as  shown  in  Fig.  7. 

Two  types  of  initial  beam  shapes  will  be  studied:  (1)  a  solid  centered 
gaussian  beam  and  (2)  a  hollow-centered  donut  beam.  We  wish  to  learn  how 
well  the  ring  fiber  will  confine  these  two  types  of  beams.  The  field 
expression  for  a  solid  gaussian  beam  takes  the  form 

-  -  (L)2 
2  V 

u  *  e 
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while  the  hollow-centered  donut  beam  takes  the  form 


u 


.  -  §  i©2 

e  2  a 


1] 


where  a  and  a  are  given  constants.  Results  of  our  computation  are  shown 
in  Figs.  8  and  9.  By  following  the  evolution  of  the  beam  intensities,  one 
may  determine  how  well  the  ring  fiber  is  guiding  the  beam.  It  can  be  seen 
from  these  figures  that  the  solid  beams  appear  to  be  better  confined  than 
the  hollow  beams,  although  the  spreading  of  the  solid  beam  energy  is  quite 
noticeable.  It  also  appears  that  simple  insertion  of  beam  energy  in  the 
high  index  region  of  the  ring  fiber  does  not  insure  good  guidance  of  the 
beam  energy.  One  may  conclude  from  this  preliminary  study  that  neither 
solid  Gaussian  beams  nor  hollow-centered  Gaussian  beams  correspond  to 
the  mode  energy  distribution  of  a  single-mode  in  a  ring  fiber.  One  should 
first  perform  the  classical  modal  analysis  to  obtain  the  mode  pattern  of 
the  single  mode  and  then  use  this  mode  pattern  as  the  initial  beam  pattern 
for  propagation  down  the  ring  fiber.  It  is  believed  that  the  use  of  ring- 
index  fiber  as  large  core  single-mode  fiber  definitely  possesses  merit  and 
should  be  studied  further.  What  we  have  demonstrated  with  our  present 
study  is  that  our  program  is  capable  of  handling  this  type  of  fibers. 

(c)  Fiber  Couplers 

One  of  the  simplest  type  of  light  couplers  is  the  fiber  coupler. 

By  placing  two  or  more  fibers  in  close  proximity  of  each  other  light 
energy  may  transform  from  one  to  the  other  through  the  coupling  effect. 
This  coupling  process  is  rather  involved.  The  well-known  coupled  mode 
theory  may  be  adequate  for  simple,  single-mode  structures  such  as  slabs 


with  reasonable  separations.  But,  when  multi-mode  complex  structures  such 


as  the  fiber  couplers  are  involved,  the  coupled  mode  theory  becomes  grossly 
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inadequate.*  On  the  other  hand,  our  FFT-scalar  wave  approach  is  uniquely 
qualified  to  deal  with  this  fiber  coupler  problem.  This  is  because  this 
technique  provides  the  evolution  of  beam  field  as  it  propagates  down  a 
complex  multimode  inhomogeneous  fiber  structure.  Four  types  of  fiber 
couplers  have  been  studied: 

Case  1  Coupling  between  two  equal  parabolic  index  fibers. 


Two  graded-index  fibers  are  fused  together  longitudinally  with 
separation  d  between  their  centers.  The  index  profile  for  each 
fiber  is  given  by 


where  n^,  6,  and  a  are  given  constants,  and  1  or  2  refers, 
respectively,  to  //I  or  II 2  fiber.  Typical  values  for  a  Corning 
or  ITT  graded  index  fibers  are  used: 

nQ  =  1.48 
6  =  0.02 
a  =  25pm 


*Recent  advances  by  L.  Eyges  and  P.  Oianino  of  RADC  using  the  extended 
boundary  condition  technique  have  shown  that  single  mode  couplers  involving 
arbitrarily  shaped  uniform  core  guides  can  be  successfully  and  accurately 
treated. 


Various  separation  d  were  used.  A  pmssian  beam  represented  by 


u(x,y)  =  uQ  exp  |  J-  (x  +  ^>2  -  y2J  /w? 

where  u(x,y)  is  the  scalar  wave  function  of  the  beam,  and  u^, 
w  are  given  constants,  is  incident  on  one  of  the  fibers.  Results 
have  been  obtained  for 

w  =  2.5u^v  5um,  lOum 
d  =  8pm,  12iin,  16pm,  20]m. 

The  evolution  of  the  beam  along  this  structure  is  shown  in 
Figs.  10-11  Displayed  in  Figs.  12  is  the  %-power  in  one  fiber 
as  a  function  of  longitudinal  distance  for  various  separations 
and  initial  beam  sizes.  For  the  situations  considered  above, 
many  modes  are  excited.  The  coupling  process  is  very  involved 
as  displayed  in  Fig.  10-11  It  still  appears  that  back  and 
forth  power  exchange  among  the  guides  prevails.  The 
complexity  of  the  power  exchange  phenomenon  for  the  multimode 
coupler  re-emphasizes  the  importance  of  obtaining  design  data 
through  analysis  before  the  actual  construction  of  fiber  coupler. 

Case  2  Coupling  between  two  equal  step-index  fibers. 

This  coupler  is  identical  to  the  previous  one  except  step-index 
fibers  were  used.  We  shall  approximate  the  index  profile  of  a 
step  index  fiber  by  the  following  expression: 
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FIG  12  (a) 


with  m  =  4 .  Again  typical  values  for  a  Corning  or  ITT  graded 
index  fiber  are  used;  i.e.,  n^  =  1.48,  6  =  0.02,  a  =  25pm. 
Results  have  been  obtained  for 


w  =  2 . 5pm ,  5  pm 
d  =  8pm, 12pm,  16pm 

where  d  is  the  separation  distance  and  w  is  the  beam  waist  radius. 

Specific  results  are  given  in  Figures  13.  Displayed  in 

Figs.  14  is  the  %  power  in  one  fiber  as  a  function  of  longitudinal 
distance  for  various  separations. 

Case  3  Coupling  between  two  unequal  fibers. 

It  is  of  interest  to  learn,  when  two  unequal  size  fibers  are  placed  side  by  side, 
whether  transfer  of  power  would  occur  for  the  multimode  case.  The  following 
fibers  were  used: 


=  1.48  m  =  1,  4 

«  =  0.02 


a^  =  25pm 
“  12.5pm 

Separation  distance,  d  *=  8^n , 12:® 
Initial  Beam  Radius,  w  ■  5wn,  10pm 
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For  the  above  chosen  parameters,  one  may  observe  from  Figs.  15-16 
that  only  nomina1  coupling  occurs  between  two  unequal  size  fibers. 

In  other  words  these  structures  are  not  efficient  couplers. 

Case  4  Coupling  between  more  than  two  fibers. 

As  illustrations,  we  have  considered  two  types  of  structures:  Three  equal 
size  fibers  located  equal-distance  from  each  other  and  three  unequal  size  fibers 
arranged  in  a  triangle  shape.  One  fiber  is  initially  illuminated,  we  wish  to 
learn  the  power  exciiange  characteristics  in  these  two  couplers.  The  following 


a^  =  25pm 

a^  =  25pm,  20pm 

a^  =  25pm,  15pm 


Separation  distances,  =  d^  =  d^  *  12pm 
Initial  beam  radius,  w  =  5um 

It  can  be  seen  from  Fig. 17  that  when  one  of  the  fiber  of  3 
identical  fibers,  separated  an  equal  distance  from  each  other, 
is  illuminated  with  a  gaussian  beam,  power  exchange  takes  place 
between  two  unilluminated  fibers  with  the  one  illuminated 
fiber  in  a  synchronous  manner.  In  other  words,  each  of  the 
originally  unilluminated  fibers  is  receiving  identical  amount 
of  power  transfer  from  the  illuminated  one,  as  expected.  On 
the  other  hand,  the  power  exchange  phenomenon  for  3  non-identical 
fibers  coupler  is  much  more  complicated  as  seen  from  Fig. IF 
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Larger  power  transfer  tends  to  take  place  among  fibers  with  s :>>•*'■>*• 

Fiber  coupler  is  one  of  the  crucial  components  in  any  fiber  optics 
system.  As  such,  one  must  understand  the  detailed  wave  interaction 
phenomenon  in  this  type  of  coupler  so  that  correct  design  can  be  made. 

We  have  demonstrated  the  capability  of  our  technique  in  dealing  with  this 
type  of  fiber  structures.  Systematic  studies  of  different  variaties  of 
fiber  couplers  may  now  be  undertaken. 

(d)  Reflection  Coefficient  Calculations. 

We  have  implemented  the  heuristic  approach  discussed  earlier  in 
the  computer  program  to  yield  reflection  coefficients  for  waves  propagating 
in  the  various  multimode  fiber  structures.  Although  it  is  difficult  to 
justify  the  accuracy  of  the  absolute  values  for  the  reflection  coefficients 
obtained  according  to  this  algorithm,  nevertheless,  we  feel  that  their 
relative  values  for  different  structures  can  be  believed.  This  is  because 
our  heuristic  approach  took  into  consideration  the  fundamental  characteristic 
of  wave  reflection:  i.e.,  reflection  occurs  when  discontinuity  of  the 
propagation  medium  or  structure  is  experienced  by  the  wave.  The  larger  is 
the  discontinuity,  the  larger  will  be  the  reflected  energy. 

Results  for  sample  calculations  for  the  reflection  coefficients  for 
various  fiber  structures  are  shown  in  Figs. 19-20 


IV.  Conclusions  and  Recommendations 

Support  of  this  program  has  enabled  the  contractor  to  develop  and 
perfect  a  computer  program  based  on  the  scalar  wave  -  FFT  algoriti  n  to 
study  the  propagation  characteristics  of  guided  waves  in  several  important, 
practical  fiber  structures  such  as  fibers  with  general  index  profiles 

(step  index,  parabolic  index,  ring  index,  etc),  multi-channel  fiber  couplers, 
and  fiber  horns  or  tapers.  These  fiber  structures  may  be  made  with  commercially 
available  fibers  whose  index  variation  may  be  as  large  as  1  -  2%. 

We  have  implemented  the  adaptive  coordinates  and  lossy  outer  mesh 
boundary  schemes  in  our  computer  program.  However,  for  most  practical 
situations  of  interest  in  which  the  fiber  core  radius  is  about  25ym,  the 
cladding  index  is  about  1-2%  less  than  the  core  index,  the  spot  size  of 
the  beam  is  less  than  20um  and  the  free-space  wavelength  of  the  beam  is 
larger  than  0.6ym, it  is  not  necessary  to  implement  the  adaptive  coordinates 
and  lossy  outer  mesh  boundary  schemes.  We  also  discovered  that  steep  index 
gradient  is  not  a  hindrance  for  the  program  to  produce  accurate  results  as 
long  as  the  scalar  wave  appraoch  is  still  justifiable. 

It  is  not  unreasonable  to  ask  the  following  question: 

"Now  that  we  have  completed  a  beautiful  program  capable  of 
producing  propagation  results  for  a  variety  of  practical 
fiber  structures,  what  can  we  do  with  it?" 

The  answer  is  "May  be  a  lot!"  Listed  below  are  only  a  few  of  the  problems 

that  we  can  solve  with  this  program: 

(1)  Any  single-mode  or  multimode  weakly  guiding  fiber  with  arbitrary 
refractive  index  profile. 

Our  program  provides  the  means  to  obtain  the  propagation  characteristics 
of  guided  waves  supported  by  this  structure.  The  core  of  the  struc¬ 
ture  may  be  circular,  elliptical,  rectangular,  triangular  or  dumb-bell 
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shape  with  general  index  profiler 

(2)  Any  fiber  couplers  composed  of  parallel  strands  of  two  or  more 
of  the  above  fibers.  This  is  the  only  program  which  can  provide 
the  detailed  coupling  characteristics  of  this  type  of  structure. 
Prior  knowledge  of  coupling  characteristics  of  a  coupler  is  the 
key  to  successfully  design  and  construct  fiber  couplers. 

(3)  Any  transition  elements  derived  from  the  above  fibers. 

Transition  elements  such  as  tapers,  horns,  or  mode  converters  or 
branches  can  all  be  analyzed  by  our  program. 

Recommended  Future  Research 

In  addition  to  the  important  practical  problems  mentioned  above  that 
can  be  solved  by  our  approach,  it  is  worthwhile  to  look  into  the  future  and 
seek  out  problems  of  potential  importance  and  interest.  For  example: 

(1)  Nonlinear  Fiber. 

Very  high  intensity  is  achievable  in  fibers.  One  may  wish  to  learn 
the  propagation  characteristics  of  waves  in  a  fiber  in  which  the 
induced  nonlinearity  of  the  medium  plays  a  significant  role.  This 
problem  may  be  solved  by  the  scalar  wave  -  FFT  approach. 

(2)  Mode  Conversion  in  a  Fiber  Due  to  its  statistically  Varying  Medium. 

This  problem  may  also  be  approached  from  the  scalar  wave  -  FFT 

point  of  view.  It  is  known  that  the  presence  of  frozen-in 

statistically  varying  medium  contributes  to  the  mode  conversion 

9 

phenomenon  in  a  fiber.  A  systematic  study  of  this  problem  will 
reveal  the  severity  of  this  effect  in  changing  the  dispersion 
characteristics  of  beam  in  this  fiber. 

(3)  large-size  Single-Mode  Fiber 

The  advantages  of  having  large-size  single-mode  fiber  are  well- 
known.  Ring-index  fiber  as  proposed  by  L  Eyges  of  RADC  appears 
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to  be  a  promising  one.  Other  type  of  fiber  whose  index  variation 
may  be  radially  unsymmetrical,  such  as,  a  layered  fiber  as  shown 
in  Fig.  21  may  also  be  promising.  Research  should  be  encouraged 
on  this  type  of  fibers. 

(4)  Polarization  Preserving  Fibers 

One  of  the  main  features  of  a  weakly  guiding  fiber  which  can  be 
analyzed  by  the  scalar  wave  approach  is  that  the  guiding  structure 
is  polarization  insensitive.  However,  for  several  important 
application  areas,  the  polarization  preserving  characteristic 
of  fiber  is  essential.  Stress  induced  birefringent  fiber  or 
deformed  core  fiber  may  satisfy  the  needed  requirement.  However, 
to  achieve  the  targeted  isolation  for  the  two  orthogonal  dominant 
modes,  the  required  stress  is  excessively  high  for  stress  induced 
birefringent  fiber  and  the  loss  is  excessive  for  deformed-core 
fiber.  We  propose  the  use  of  layered  dielectrics  as  shown  in 
Fig.  21  to  produce  an  equivalent  birefringent  effect  to  enable 
proper  isolation  for  the  two  orthogonal  dominant  modest®  Initial 
indication  is  very  promising.  It  would  therefore  be  very  worthwhile 
to  pursue  this  research. 


1 
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FIG.  21  Polarization  Preserving  Layered  Fiber 
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SAMPLE  PROGRAM  LISTINGS  FOR  THE  CASE  OF  BEAM 
PROPAGATION  IN  A  STEP  INDEX  FIBER 


CORE  SIZE  (DIAMETER)  =  50ym 
CORE  INDEX  =  1.48 
CLADDING  INDEX  =1.46 

INITIAL  GAUSSIAN  BEAM  SPOT  SIZE  (DIAMETER)  =  40pm 
WAVELENGTH  =  0.8pm 

(STEP  INDEX  CASE,  m  =  6) 
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SRFT-SCALE  CHARACTERS  ANO  RANGES 


0.5976  OAE- 02 

0.1  053  78E  0  0 

• 

0.1053  78E400 

0.2  04  78 1 FF 00 

— 

0.204781  E+-  00 

0.304  1  83  E  0  0 

♦ 

0.3041  83F  f 00 

0. 403585F+00 

* 

0.403586F*  00 

0.5029  88F  0  0 

ff 

0.5029  38E+00 

0 .602390F+00 

» 

0 .6023  90S  00 

0.  701  793E#-  00 

• 

0. 701 793E400 

0.801 l 95E+00 

■ 

0.8011 95F400 

0.900598E*  00 

■ 

0.9005  93E4  00 

0.1 OOOOOEFO  1 

THIS  IS 

STEP  1 

SUML  = 

0. 4859293EA00  SUNR 

=  0.51 40707EA00 

TEAM  WAIST  IN  X  ANO  V  DIRECT 

IONS  IS  19.50358 
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THIS  IS  STEP  2 
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5 RET  — SC  ALE  CHARACTERS  AND  RANGES 


0.323949E-02 
0.21  t  I  44E+00 
0.41 90  48E  *  00 
0.626953E+00 
0. B34857E*  00 
O.  1042  76r  +  0l 
0.1 250  6  7F  ♦  0  1 
0.14SB57E*  01 
0. 166647E*01 
0. 18743BE401 


0*2111 44E*  0  0 
0.41 9048E+0  0 
0.626953E*  00 
0.8348S7E*00 
0.1042  76E*0  1 
0. 1 2S066E*  0  1 
0*1 4SR57F  *01 
0.  1  6664  7E  *01 
0.1 87438E+0 l 
0*20  82  2  8E  *01 


THIS  IS  STEP  3 


SUNL  *  0.48171 33E *00  SUMR  *  0 . 51 82866E *0 0 

SEAN  MA 1ST  (N  X  AND  T  DIRECTIONS  IS  18.07376  18.07498  MICRONS 
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GRET-SCALE  CHARACTERS  AND  RANGES 


0.7596  99E-05 
0  •  842067F.  ♦  00 
0.  1684  1?E«-Cl 
0.2526  18E4-01 
0  •  3368  24E  01 

0.4210306*01 
0.5052  366*  01 
0 .5894  4  26*  0  t 
0.6  736  4  86*  01 
0 • 7578  54  E*  0 1 


THIS  IS  STEP  4 


0.S42067E*00 
0.1 684 l 3F ♦ 0 l 
0.2  526  I  8E  *01 
0.336824E*0 I 
0.  42  1  030E*01 
0.5052  36F  *01 
0 .5894  426+-0  I 
0.67  36  48F *01 
0.  75  7054E  *  0  I 
0 .8420606*0  I 


SUNC  »  0.4  754  75  06*00  SUMR  *  0 . 52  4S2  4  9E +0  0 

BEAM  WAIST  IN  X  AND  V  O  IRECT IORS  IS  17.64449 


17.04572  MICRONS 
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jaey-$cXi_E  Characters  and  ranges 
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0.306083F*  Cl 
0.4591  12F*0  l 
0.6121  42F4-  0 1 
0. 7651  7  IF  ♦■01 
0.918200F+01 
0.1071  23E*02 
0. 122426E* 02 
0.1 37  7  29E ♦ 02 


0. 153054E*01 
0.306083E*01 
0.4591 12E*01 
0.6121 4?F  *01 
0.7651 71Ff01 
0.91 8?00F*0  l 
0.  10  71 2  3E*  02 
0. 1 224  26E ♦ 02 
O. 1 37729E+02 
0. 153032E+02 


IS  STEP  5 

=  0  •  4 803946E  +  00  SU4R  =  0 . 51  96053E  *■  00 

WAIST  IN  X  ANO  V  DIRECTIONS  IS  17.43616  1 7. 43765  MI  CRONS 
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SREY--SC  ALE  CHARACTERS  AND  RANGES 

0.34857BE- C5 
.  0.5467  19E*00 

0. 109343F* 01 

•  0  a  1 640  1  5E *  0  l 

•  0.218bS6E+01 

0  0.273359F*01 

■  0.328029E4-01 

S  0  a  3827  OIF*  01 

•  0.437373E*0I 

n  0  a  4920  4 4F  *  C 1 


0  a  5  4  6  7  1  9E  *00 
0. I09343F+0 1 
Oa  1540 I5E*0 1 
0  a  2 1  8ft  8ftF  *01 

0.273358E+0 1 
Oa  32  802  9E«-  0 l 
0.3«270lE*0l 
Oa  4  3  7  3  72 F  *  0  l 
0.492044E40 1 
Oa  54671 6E*  01 


THIS  IS  STEP  6 


SURL  *  0,4  704 280 E*  00  SURR  *  0.52957 19E* 00 

3 EAR  WAIST  IN  X  AND  Y  DIRECTIONS  IS  17.47115  17.47259  MICRONS 
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SRE  Y  — SC  ALE  CHARACTERS  AND  RANGES 


* 

* 

0 

K 


THIS  IS  STEP 


0.341S98E-C4 
0.  146  743F *01 
0.293483E*  01 
0.44022  ir*0  1 
0. 586963E*  0 1 
0.7337  03E  *01 
0.880443E+01 
0. 1027 1 8F ♦ 02 
0.  I  1 7392E*  02 
0.1 320  66E  *  02 


0. 146743E+0  1 
0.293483E+0 l 
0.440223E*  01 
0.586963F*01 
0.  7  13703E*0 1 
0.  88044 3 E  *■  0  1 
0 . 1  02  7 1 8E*  0  2 
0.  I  1 7392E*  02 
0.  1 32066E  +  02 
0.  146740E  *02 


S URL  *  0.4  73628 1 F* 00  SCMR  =  0 . 5263 7 1 8E *0 O 

3EAR  WAIST  IN  X  AND  Y  DIRECTIONS  IS  17.7574S 


17. 75934 
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0.3495 78E-03 

0.996981 E*  0  0 

. 

0.9969  81E*  00 

0. 199361F*01 

~ 

0.  199361E*01 

0.2  990  24  C  0  1 

* 

0.299024F*  01 

0.398687E+0 1 

* 

0.398688F*  Cl 

0.  498351  E 01 

0 

0  .498  3  5 1  E  0  1 

0.5980 14F*0  l 

« 

0.5980  I4F>  01 

0. 6976  77E«-0  1 

e 

0 .6976  77E*  0  1 

0. 797340E*0  1 

• 

0.797340F*  01 

0. 897003E+0I 

■ 

0.8970  03F  *0 1 

0 ■ 996667E* 0  I 

THIS  IS 

STEP  8 

SUML  * 

0.4  80O833E*  00  SUMR  = 

0 .51 991 67E400 

3EAM  WAIST  IN  X  ANO  Y  DIRECTIONS  IS  18.20370  18.20607  MICRONS 
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0.753595E-04 
0  •  3  798  3  OE  *0  0 
0.7595S4E4-00 
O. 1 l 39  34E ♦ 0  1 
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0. 3 79R30E4  0 0 
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r  HI S  IS  STEP  10 
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0. 10321  3E-  C3 

0.410842E+00 

. 

0.*l  08*2E«-00 

0  •  82  1  58 1  £♦■  0  0 

- 

0.821582E#-  CO 

0. 1 23232E*  01 

♦ 

0 • 1 232  32  F ♦01 

0.  16*3 06E 0  l 

* 

0  •  116*3  06  E  +  0  1 

0 . 20  S3  80E*  0  l 

A 

0.2053  80Ef  01 

0. 246454E*0l 

« 

0.2*64S*F>01 

0.28752  BE  *  0 1 

• 

0. 28752 8E f  01 

0.328  60  2  E  ♦ Ol 

• 

0.3286  02E+-01 

0.  J696  76E«-0l 

• 

0.3096  76E*  01 

0.*1  0749F«-0l 

THIS  IS 

STEP  11 

3  URL  a 
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SREV  — SC  ALE  CHARACTERS  AND  RANGES 

O. 3 73 785 E—  OA 

*  0.5941 57E* 00 
0.1188286*01 

*  0*1  78?  39E *01 

*  0.2376516*01 

*  0.2970636*01 

■  0.3564756*01 

*  0.4158876*01 

■  0  •  A  752  99E  *  0 1 

■  0.S3471 1 E*OI 


THIS  IS  STEP  12 


SUNL  *  0.4  75021 76*  00  SUMR  =  0 . 52 497B3E *00 


SEAM  WAIST  IN  X  AND  V  DIRECTIONS  IS  18.92630  18.92423  MICRONS 


0.594157E*00 
0. 1  I 8828E+01 
0. I  78?  39E ♦ 0 1 
0.23/6SlE*01 
0.297063E*  01 
0 . 3564  75E*  0 1 
0. 4158876*01 
0.4752996*0 l 
0.53471  1  E*  01 
0.5941236*01 
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5 REV -SC  ALE  CHARACTERS  AND  RANGES 

0.1262S7E-C4  0.898516E* 00 

.  0.8985 16E*00  0.179702E*01 

I  -  0.1797  02  E  0 1  0.269552E+0I 

I  ♦  O  .  2695  52E  0  l  0 . 3594  03E  ♦- 0  l 

!  *  0.359403E+-01  0 . 449253E* 0 l 

'  0  0.44925TE*  01  0.539I03E*01 

*  0.539103001  0.628954B«-0l 

1  •  0.5289540  01  0.7I8804EV01 

■  0  -  71  88  04001  0 .80  86 54E  +  0  I 

■  0.808654E*  Cl  0.  09  8505E *-01 

i 

1  THIS  IS  STEP  13 


5UNL  =  0.4702  22  7E  ♦  0  0  SUM*  =  0 . 5297772E  *■  0  O 


SEAR  HA  1ST  IN  X  A NO  V  DIRECTIONS  IS  18.73090  13.73331  MICRONS 
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0.703742E-C4 
0  .8480  2 l F  *  00 
0  •  t  695  97E  *  0 l 
0  •  254.1 92F  *01 
0.3391 87P ♦ 0 1 
0.4239R2E*0t 
0.5087  73f  *  01 
0.5935  73E ♦01 
0  .67836RE  *0 l 
0. 763 1 63E*  0 1 


0.848021 E*  00 
0. 1 695  9  7F  ►  0  l 
0.2543926*01 
0. 3391  87E*  01 
0.4239826*0  1 
0.50877BE*0  l 
0. 5935736*  01 
0.fc78368E*0  1 
0. 7631 63E*  0  1 
0.84  7958E*0  1 


THIS  IS  STEP  14 


SUWL  *  0.46514386*00  SUMR  =  0. 534  856 1  E*0  0 

SEAN  WAIST  IN  X  ANO  V  DIRECTIONS  IS  18.77950  18. 


76253  MICRONS 
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3REV-SCALE  CHARACTERS  AND  RANGES 


0.2441 I4E-03 

0.1 55652E«-Q  1 

« 

0. I55652F4  01 

0.31  1279E+01 

- 

0.31  1279E4-01 

0.4669  06  E  0  1 

♦ 

0.4669  06F*  0 1 

0  .F  225  3  3E*-0  1 

• 

0.622534F*  01 

0 •  77 8 1  6 1  E  ►  0  l 

0 

0.7781  6 1  E  0 1 

0.933  7  88E ♦01 

* 

0.933  7  88F  *01 

0.  1  08942E  #-02 

• 

0.  108942F4C2 

0. 124504E4-02 

• 

0.1  245  04E  02 

0. 1 40067E*  02 

• 

0.140067E*02 

0.1 55630E402 

THIS  IS 

STEP  IS 

sunl  « 

0. 4  69  0364  E4  00  SUMR  = 

0.53096 35E+00 

SEAM  WAIST  (N  X  AND  Y  DIRECTIONS  IS 


18.87981 


i a. aaa 


F/6  20/6 


1 


MICROCOPY  RESOLUTION  TEST  CHART 


GRET-SCALE  CHARACTERS  AND  RANGES 


0*2469 0SE- 03 
0*6976  66E4-00 
0*I395  08E>01 
0*2092 50F4 01 
0.276992E+01 
0.348734E+01 
0*4184  76«r*0l 
0.488? 18E+01 
0  *55 79 6 OF 40  t 
0.627702E401 


0.697666F400 
0* 1 39509E4  0  I 
0.2092S0E  +  0  1 
0*  278992F.401 
0.348734E401 
0*41 S476E+01 
0*4  882 l  8E ♦ 0 1 
0.557960E40  1 
0*  627702E  *01 
0.697444E40  1 


THIS  IS  STEP  16 

SU«L  *  0.4632549F *00  SUMR  *  0 • 536 74 50E40 0 

SEAM  WAIST  IN  X  AND  V  DIRECTIONS  IS  19.10669 
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0.91 I 394E— 03 
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0.5264 30E4’  01 
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0.  1  05?  77E4  02 
0.1 18435E402 


0.  1  316  76E 4-01 
0 .26326 1 E  +  0 1 
O. 3948 4 5E ♦ 0 1 
0.5264S0E401 
0.6580 15E40  1 
0. 789600E*  0 1 
0.921  1 85E  *01 
0.105? 77E402 
0.1184  35E+02 
O.  1  31594002 


THIS  IS  STEP  17 


SUNL  *  0 .461 8073E «-00  SIM*  *  0.S381926E4-00 

9EAM  MA  1ST  IN  X  AND  V  DIRECTIONS  IS  19.16985 


19.17290  MICRONS 
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3REY-SCALE  CHARACTERS  AND  RANCCS 

0.192583E-04 

•  0.  1595426401 
0.  31 9081 64  01 

4  0  *4  78621  F  ♦•0  1 

•  0.6381 616401 

0  0.79 770064  01 

ft  0.9572406401 

■  0.1116  786  4C2 

•  0.  1276326402 

•  0.143586F402 


THIS  IS  STEP  18 


SUMl  «  0.47309606400  SUMR  *  0.52690326*00 


36AM  VAIST  IN  X  AND  T  DIRECTIONS  IS  19.34  006  19.343  09 


0.159S42E401 
0.31 908 1 E  +  0  1 
0.4  7 86? IE *01 
0.638161F401 
0.797700E40 I 
0.95724  OF  A  0 1 
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0. 127632E402 
0.1435S6E402 
0. 1595406402 
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